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Procedure for Separating Noise Sources in Measurements of 
Turbofan Engine Core Noise 

Jeffrey Hilton Miles 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


Abstract 

The study of core noise from turbofan engines has become more important as noise 
from other sources like the fan and jet have been reduced. A multiple microphone and 
acoustic source modeling method to separate correlated and uncorrelated sources has been 
developed. The auto and cross spectrum in the frequency range below 1000 Hz is fitted 
with a noise propagation model based on a source couplet consisting of a single incoherent 
source with a single coherent source or a source triplet consisting of a single incoherent 
source with two coherent point sources. Examples are presented using data from a Pratt 
and Whitney PW4098 turbofan engine. The method works well. 


Nomenclature 


A amplitude 

B amplitude of a point source B 

B e resolution bandwidth, Hz., B e = 1/Td = r/NP = 11.71875 Hz 
c 0 speed of sound, ml sec 
E expected value 
/ frequency 

F\\ two sided Fourier transform operator 

f c upper frequency limit, f c = l/2At = r / 2, Hz. (24000 Hz.) 

G xx (f) auto power spectral density function defined for non- negative frequencies only (one-sided) 
G xy (f) cross power spectral density function defined for non- negative frequencies only (one-sided) 
j positive imaginary square root of — 1,\AM 

k wavenumber 

L y number of frequencies, / C /A/ = A/2 ( 2048 ) 

MW molecular weight of air, 0.02897. 

b ’ . ’ mol 

A amplitude of uncorrelated noise source 

NP segment length, number of data points per segment (4096) 

rid number of disjoint (independent) segments , rid = B e T to tai = 234 

n 0 number of overlapped segments/blocks 

p(t) pressure wave function 

Pi probability confidence interval, percent 

R gas constant, 8.314, 

r radial distance, r 

r sample rate, samples/sec. (48000) 

R xx {t) auto-correlation function 

R xy {r) cross-correlation function 

S xx two sided auto spectral density function 

S xy two sided cross spectral density function 

Ttotai total record length, sec. (« 20 sec.) 
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t time, sec. 

Td record length of segment 

TF air temperature degrees Fahrcnheit(80°) 

TK temperature in degrees Kelvin 

X(f) two sided Fourier transform 

X$(i, f) one sided Fourier transform of block i 

Subscripts 

o some arbitrary specific value 

x signal x 

y signal y 

Symbols 

A / frquency step, l/X^Hz. (11.718) 

At sampling interval, 1/r (1/48000), sec. 

5 dirac delta function 

7 ratio of specific heats, 1.4 

7 nn(f) magnitude squared coherence (MSC) function of noise 
7 2 y (/) magnitude squared coherence (MSC) function 
r time displacement, sec. 

T\ propagation time delay 

I. Introduction 

Understanding turbofan engine noise is a key priority of the National Aeronautics and Space Administra- 
tion (NASA). Consequently, new diagnostic procedures to identify dominant sources and changes in dominant 
sources are being developed. This paper discusses a new diagnostic procedure that separates correlated far 
field turbofan noise from uncorrelated farfield turbofan noise using a multiple microphone method and point 
source acoustic models. The auto and cross spectrum in the frequency range below 1000 Hz is fitted with 
a noise propagation model based on a source couplet consisting of a single incoherent source with a single 
coherent source or a source triplet consisting of a single incoherent source with two coherent point sources. 
Two noise signals are coherent if they can be aligned with one another so that the coherence calculated 
using a periodogram averaging method is greater than the coherence of two random signals using the same 
periodogram averaging method. As fan noise and jet noise from turbofan engines are reduced, the issue 
of core noise reduction becomes more important. While core noise may be reduced by acoustic liners, the 
issue of measuring the amount of reduction becomes increasingly important. The proposed scheme separates 
coherent noise from random jet noise and should enable the measurement of the effectiveness of core noise 
reduction liners. The scheme was developed as part of a research program to study core noise from a Pratt 

6 Whitney PW4098 turbofan engine ( Miles. 1-3 ) 

Procedures using coherence-based techniques have been developed for extracting acoustic signals buried 
in noise. The coherent output power method for noise source identification is discussed in by Bendat. 4 ~ 6 The 
application of this technique that is of interest is the use of coherent output power spectra to separate 
and identify correlated combustion noise in far field measurements of turbofan engine noise. Karchmer 7 
and Karchmer, Reshotko, and Montegani, 8 use the coherence function calculated from internal microphone 
measurements of fluctuating pressures in the combustor and far field acoustic pressures to determine the 
correlated combustion noise of a YF102 turbofan engine at far field locations by calculating the coherent 
output power spectrum. 

The three signal coherence technique was developed by Chung 9, 10 for flow noise rejection. A similar 
technique was developed and used by Krejsa. 11 The three signal coherence technique was used by Shiv- 
ashankara 12 to study core noise in a Pratt and Whitney JT9D. It was used by Hsu and Ahuja 13 to separate 
ejector internal mixing noise from far field measurements and by Stoker, Ahuja, and Hsu 14 to separate 
wind-tunnel background noise and wind noise from automobile interior noise measurements. It was used 
by Michalke, Arnold and Holste 15 to study sound in a circular duct with mean flow. The method is also 
discussed by Minami and Ahuja. 16 

The inherent coherent properties of jet and core noise have been used in several source separation proce- 
dures. The radiated field of jet noise is has a low coherence when measured at two widely separated points. 
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Parthasarathy 1 ' attributes this to the fact that the jet noise sources are in motion so that the source fre- 
quencies undergo large Doppler shifts as the noise is radiated to the far field. Core noise has a high coherence 
when measured at two widely separated points. This is attributed to the fact that the frequencies of the 
radiated core noise are preserved unchanged in the far field Parthasarathy. 17 Consequently, three signal 
coherence technique is especially useful to study turbofan core noise using widely spaced microphones. 

A method of separating jet noise and core noise using auto-correlations and cross-correlations was devel- 
oped by Parthasarathy. 17 The model presented assumes the microphones are located on the arc of a circle 
about a single source having a known position. This is a couplet source which produces a core noise sound 
radiation field that is correlated over wide microphone spacing and a jet noise sound radiation field that is 
not correlated over wide microphone spacing. Since the microphones are on the same arc, the model cross- 
correlations calculated have zero propagation time delay. In addition, a model for the jet noise directivity is 
used. The resulting system of three unknowns and three equations can be solved exactly. In applying this 
method to data from a small jet the source location is known and a correction to the amplitude correlation 
measurements is made so that the microphones are located on the arc of a circle about a source having a 
known position. In applying this method to measured data, I assume the cross-correlations were also time 
shifted to remove the propagation time delay. An extension of the method to obtain spectral information 
is also discussed using non zero values of propagation time delay. The phase angle shift due to propagation 
time delay is not discussed. Again, this seems to indicate that the propagation time delay is being set to 
zero. The method calculates the jet noise spectrum, the core noise spectrum, and the directivity ratio of 
the core noise at a particular angle as a function of frequency from measured auto and cross spectra. The 
method was applied to experimental data obtained in a small scale facility. The spectral formulation of this 
method was used by Tessson 18, 19 to study jet noise and core noise from the static test of a small gas turbine 
engine in an anechoic chamber. 

A method to achieve separation of two different correlated noise sources from far field measurements 
contaminated by extraneous jet noise using multiple microphones was developed by Minami and Ahuja. 16 
The equations discussed use measured auto-spectra and cross-spectra. The model assumes the source noise 
can be represented by a triplet consisting of a correlated noise source A, correlated noise source B, and 
an uncorrelated noise source. I assume that all the sources are at a single triplet location however in the 
formulation presented no source location or microphone location information is used. The five microphone 
method described involves solving a set of 55 equations for 55 unknowns at each frequency of interest. 
At each microphone location the spectrum of correlated noise source A, correlated noise source B, and an 
uncorrelated noise source is obtained as a function of frequency. The method was validated using analytical 
simulations. 

The method presented herein also models the auto- and cross- spectrum measurements made with multiple 
microphones. The method uses a point source propagation model so the microphones need not be on an arc. 
As a consequence, the directivity is the same in each direction for each source. The method assumes the 
engine noise source can be modeled with a source doublet or triplet. One part source doublet and triplet 
consists of an incoherent source at the location taken as x = 0.0 and y = 0.0. The source doublet has an 
additional coherent source at location x = xa with y = 0.0. The source triplet has two additional coherent 
sources at location x = xa with y = 0.0 and x = xb with y = 0.0. As part of the procedure a comparison of 
results obtained with the use of a source couplet with the results obtained using a source triplet is made. In 
addition, a source position is assumed to have a y coordinate of zero and the x coordinate along the engine 
axis is found as a function of frequency. This is done since noise can be radiated from the engine case to 
the far field as well as from the nozzle. In order to solve the resulting acoustic equations for a small set of 
parameters, a solution method that provides an optimum solution in a least square sense without derivatives 
was used. The method is set up to do a least squares curve fit to match the auto-spectra and the cross spectra 
magnitude and phase measurements between all microphones. The method assumes the existence of one or 
more point sources that can produce the same measurements. The problem then becomes one of finding the 
appropriate point sources. The current strategy is to have the curve fit adjust the amplitude and x location 
of the sources that determine the cross spectrum at each microphone and the amplitude of a random noise 
source. For the single coherent source case, three parameters are found at each frequency. When a model 
with two coherent sources is used, one has five parameters to determine at each frequency. Examples are 
presented calculated with a four microphone array using data from a Pratt and Whitney PW4098 turbofan 
engine. 
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PW4098 Engine at C-ll Stand for EVNRC Phase 2 Tests 



Figure 1. Pratt &: Whitney test stand Cll, West Palm Beach Florida, with PW4098 engine and attached 
acoustic inflow control device also with and without aft acoustic barrier walls for EVNRC Phase 2 tests. 
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II. Analysis Method 


The core noise is assumed to be propagating in 
acoustic modes in the turbofan engine. The co- 
herent acoustic energy leaves the nozzle and travels 
through the turbofan engine shear layer to a ground 
microphone. Additional acoustic energy from the 
jet and from random scattering reaches the same 
microphone. In this test program four microphones 
at 150 feet and angular position of 100°, 110°, 120° 
and 130° measured from the inlet were used. Con- 
sequently, the measurements available at each test 
condition are four sound pressure spectrums and six 
sound pressure cross-spectrums each consisting of a 
magnitude and a phase angle. These sixteen mea- 
surements are available as a function of frequency. 

Four real variable and six complex variable acoustic 
equations relate the measurements and the model 
parameters. The angular separation of at least 10 
degrees means that the jet noise measured at any 
microphone can be assumed to be incoherent with 
the jet noise at any another microphone. 

The basic procedure is to assume that some 
combination of coherent sources and an incoherent 
source will produce the measured auto and cross 
spectra. The coherent source is assumed to be on 
the turbofan engine centerline as some x position to 
be determined with y = 0. The incoherent source 
is assumed to be at x = 0 and y = 0. Using the 
results in appendix A on page 11 where the one- 
sided spectrum functions for acoustic signals from a 
point source are derived, the acoustic model equa- 
tions shown in appendix B through E are defined. 

The following models were considered. 

• Two parameter model. Coherent source of magnitude A at (x = x a , y = 0) (Appendix B on page 12 ). 

• Three parameter model. Coherent source of magnitude A at (x = x a ,y = 0) (Appendix C on page 13)- 
with an incoherent source of magnitude N at (x = 0, y = 0) 

• Four parameter model. Coherent source of magnitude A at (x = x a , y = 0) and a Coherent source of 
magnitude B at (a; = Xb, y = 0) (Appendix D on page 13 ). 

• Five parameter model. Coherent source of magnitude A at (x = x a , y = 0) and Coherent source 
of magnitude B at (x = Xb,y = 0) with an incoherent source of magnitude N at (x = 0,y = 0) 
(Appendix E on page 13) 

Consequently, the number of equations greatly exceeds the number of model parameters used in the study 
discussed herein. Note that no radiation pattern is assumed for these models. These are point sources that 
radiate sound equally in all directions. 

The measured auto-spectrum and cross-spectrum required as input to the acoustic model equations is 
experimentally determined and subject of nominal experimental error and statistical uncertainties. Also, 
the measured auto-spectrum includes random noise in addition to coherent signals from propagating waves. 
Consequently, a solution method that provides an optimum solution in a least squares sense without deriva- 
tives was used. Algorithms for minimization without derivatives are discussed by Brent. 20 The search 
technique used in this study is described by Powell 21 and Fortran computer code for this algorithm is given 
by Shapiro 22 and Kuester. 23 The code used was a modified version of the one in Shapiro 22 which was updated 
to be in a FORTRAN 90 style. The cost function used is written in terms of sound pressure level and phase 
angle. 
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Figure 2. Acoustic arena and microphone array at 
Pratt & Whitney test stand Cll, West Palm Beach 
Florida for EVNRC Phase 2 tests. 


III. Experiment 


To demonstrate the usefulness of the procedure for separating correlated and uncorrelated noise sources 
measurements made in the far field of a Pratt & Whitney PW4098 turbofan engine will be used. The 
measurements were made in a study of aircraft turbofan engine core noise conducted as part of the NASA 
Engine Validation of Noise Reduction Concepts (EVNRC) Program. 

The spectral estimate parameters are shown in table 1 on the next page The signal processing algorithms 
used were written in Fortran. They are based on subprogram modules developed by Stearns and David. 24 In 
the calculations the segments were overlapped by 50 percent. Fig. 1 on page 4 shows the test stand. Fig. 2 on 
the page before shows the angular placement of the far field microphones on a 150 foot radius. This analysis 
uses the microphones at 100°, 110°, 120° and 130° measured from the inlet. The coordinate system used 
herein has the x axis along the engine centerline increasing to the right. The y axis is perpendicular and 
increases toward the top of the page. The center (x = 0, y = 0) is at the engine center. 

IV. Results 


A. Cross-Spectra Validity 

The signals from each pair of microphones used to calculate the cross-spectra only produce valid measure- 
ments if the coherence is greater than some threshold coherence. The coherence function is given by 


7^ e 2 (/) — 


l G ei02(/)l" 


GW/)GW/) 

The measured coherence calculated using segments overlapped by 50 % is given by 


(1) 


70X02 


(/) = 




( 2 ) 


Er=ii^ 1 (b/)i 2 E”=ii^(b/)i 2 

In Miles 1 comparisons were made of a coherence threshold based on aligned and unaligned coherence and one 
based an analytical coherence threshold using computer simulation. Results were obtained from computer 
simulation that show good agreement with the theoretical estimate of the analytical coherence threshold 


7nn = 1 — (1 — P) 1 /^- 1 ) (3) 

where we use herein P = 0.95 and instead of the number of independent segments/blocks rid we take 
nd = n 0 where n 0 is the number of overlapped segments/blocks. The coherence threshold y 2 ra is discussed 
by Carter, 25 ’ 26 Halliday et. al. 27 (page 247), and Brillinger 28 (page 317). The coherence threshold 7 2 n has 
a value which is greater than 95 % of the values of the coherence of two independent time series calculated 
using rid disjoint data segments/blocks. The coherence threshold used herein is calculated using the number 
of overlapped segments, n 0 . As part of the set of figures for each case used to show the results, the coherence 
and the threshold coherence, 7 2 ra , will also be shown. The data used herein had coherence values above the 
threshold coherence, 7 2 „. 


B. Cross-Spectra Phase Angle Sampling Errors 

In Bendat 6 and in Piersol 29 the random error in the phase estimates due to statistical sampling is given in 
terms of the standard deviation of the estimated phase angle, 6 12, by 




' 012 (f) 

1 

~ sin -1 < 


t 1 -7 1 2 (/)] 

|7i2|-\/2?E 


1/2 ■ 


( 4 ) 


where <7 O 12 (/) is measured in radians and as used herein n 0 is selected to be the number of overlapped 
segments or blocks used in the spectral calculations. For the special case where the term in curly brackets 
is small Eq. 4 becomes 
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where for the unknown coherence 7 f 2 (/) the estimated coherence 7 g g (f) from Eq. 2 is used. A plot of the 
standard deviation of the phase angle in degrees verses coherence is shown in 3 on page 10. 

When the coherence is greater than 0.15, Fig. 3 shows the standard deviation should be less than 5 
degrees. Straight line fits was made to the phase verses frequency data. The standard deviation of the error 
was greater than on would calculate from the coherence i.e around 10 degrees when one calculated 5 degrees. 
Consequently, the phase angle measurements might be showing propagation effects due to wind shear or 
temperature gradients or a change in source location with frequency. In addition, the error verses frequency 
plots seemed to be in phase for each pair. As a consequence, instead of a using a fixed source location in the 
model the y co-ordinate of a source is fixed at zero and the x co-ordinate is free to vary. 


C. Point source models 

The method discussed has been tried out over a range of operating conditions. For each case all four models 
were used. Four typical cases presented herein shown in Fig. 4 on page 14, Fig. 5 on page 15, Fig. 6 on page 16, 
and Fig. 7 on page 17. are for N1 Corr. values of 582 rpm, 1622 rpm, 1900 rpm and 2743 rpm. Only one 
acoustic model result is shown for each case. The selected case uses the fewest number of parameters to 
achieve a good result. For each case all four auto-spectra and six cross-spectra are calculated. However, only 
the following items are shown herein: 

• The measured and calculated auto-spectrum at 100 degrees. 

• The measured and calculated cross-spectrum magnitude between the 100 degree and 120 degree mi- 
crophones. 

• The measured and calculated cross-spectrum phase angle between the 100 degree and 120 degree 
microphones. 

• The measured coherence and analytic threshold coherence between the 100 degree and 120 degree 
microphones. 

• The sound pressure amplitude curve fit parameter or parameters. 

• The correlated source location or locations. 

In the modeling scheme used herein, for 
each case the jet noise is assumed to be 
uncorrelated between the microphones 
at all frequencies. Furthermore, it is as- 
sumed the basic radiation pattern of the 
correlated noise as specified by the cross- 
spectrum phase angle could be repre- 
sented by a single point coherent source 
at some x location. In all cases, a better 
result is obtained by using an acoustic 
model which assumes a non correlated 
noise source is present and letting the 
computer solution procedure determine 
if its value is significant. In some cases 
a second coherent point source gives a 
slightly better result. Only the coher- 
ent sources create the cross spectra. The 
auto spectra are due to the coherent and 
incoherent sources. 

The general formulation is given in 
Appendix A on page 11. The general 

models are shown in Appendix B on page 12 through E on page 13. 


Table 1. Spectral Estimate Parameters 


Parameter 

value 

Segment length i.e. Data points per segment, NP 

4096 

Sample rate, r, samples/second 

48,000 

Segment length, Td = NP/r, seconds 

0.08533 

Sampling interval, At = 1/r , seconds 

2.0833 x 10 -5 

Frequency step, A / = 1 /Td, Hz 

11.718 

Upper frequency limit, f c = 1/2A t = r/2,Hz 

24000 

number of frequencies, Ly = / C /A/ = NP/ 2 

2048 

Time delay, r = 6323/48000, seconds 

0.1317 

number of independent samples 

234 

overlap 

0.50 

Sample length, sec. 

20 
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1. N1 Corr = 582 rpm 


At the N1 Corr 582 rpm operating condition the core noise is dominated by a low frequency combustion 
tone at 327 Hz as shown in Fig. 4a. The core noise radiation pattern is well represented by a single point 
source as shown in Fig. 4b and Fig. 4c. The results shown are for the three parameter couplet source 
model (A,N,x a ) described in Appendix C on page 13. As shown in Fig. 4a, at a few low frequency points 
the computer solution included a uncorrelated noise source. The coherence shown in Fig. 4d is above the 
coherence threshold = 0.00639431. Consequently, the cross spectrum phase angles look good. At most 
points the uncorrelated noise level was to low to plot as shown in Fig. 4e. As shown in Fig. 4f, in order 
to match the measured phase angle the computer model puts the source at more than x a = 50. feet for 
frequencies less than 200 Hz. and near Xb = 20 feet for frequencies greater than 200 Hz. 

2. N1 Corr = 1622 rpm 

At the N1 Corr 1622 rpm operating condition, results are shown for the five parameter triplet source model 
(A,B,N,x a ,Xb) described in Appendix E on page 13. As shown in Fig. 5a this model does a good job 
below 300 Hz. in representing the 100 degree microphone auto-spectrum and a satisfactory job at higher 
frequencies. At higher angles the curve fit at these frequencies is better. The core noise radiation pattern is 
well represented by a two point sources as shown in Fig. 5b and Fig. 5c. Below 300 Hz. this model does a 
fine job. The coherence shown in Fig. 5d is above the coherence threshold = 0.00639431 below 800Hz. 
Below 400 Hz. the coherence is grater than 0.15. Consequently, the cross spectrum phase angles look good 
below 400 Hz. The variation of the two-point model correlated source amplitude shown in Fig. 5e and the 
source separation results shown in Fig. 5a have the same frequency dependence that the combustion modes 
have in that they seem to be in 200 Hz bands. This suggests the correlated noise is due to the combustor 
source. As shown in Fig. 5f, the source locations of the two point sources have an interesting symmetry. 
Source B is near xb = 0.0 and may represent sound radiation from the engine case. Source A is near xa = 25 
and may represent sound radiation from the nozzle. 

3. N1 Corr = 1900 rpm 

Again, at the N1 Corr 1900 rpm operating condition, results are shown for the five parameter triplet source 
model ( A 7 B,N : x a ,Xb ) described in Appendix E on page 13. As shown in Fig. 6a this model does a good job 
below 200 Hz. in representing the 100 degree microphone auto-spectrum and a satisfactory job at higher 
frequencies. At higher angles the curve fit at these frequencies is better. The core noise radiation pattern is 
well represented by a two point sources as shown in Fig. 6b and Fig. 6c. The coherence shown in Fig. 6d is 
above the coherence threshold = 0.00639431 below 800Hz. At most points below 800 Hz the coherence 
is grater than 0.10. Consequently, the cross spectrum phase angles look good below 800 Hz. The variation of 
the two-point model correlated source amplitude shown in Fig. 6e and the source separation results shown in 
Fig. 6a again have the same frequency dependence that the combustion modes have in that they seem to be 
in 200 Hz bands. This again suggests the correlated noise is due to the combustor source. Again, as shown 
in Fig. 6f, the source locations of the two point sources have and interesting symmetry. Source B is again 
near xb = 0.0 and may represent sound radiation from the engine case. Source A is again near xa = 25 and 
may represent sound radiation from the nozzle. 

4- N1 Corr = 2743 rpm 

The highest power setting examined herein is the N1 Corr. 2743 rpm operating condition. The auto-spectrum 
at 100 degrees is shown in Fig. 7a. Below 100 Hz. the uncorrelated noise source is not needed to represent 
the measurements. The core noise radiation pattern is well represented again by a single point source as 
shown in Fig. 7b and Fig. 7c. The results shown are for the three parameter couplet source model( A,N,x a ) 
described in Appendix C on page 13. As shown in Fig. 7a, at frequencies above 100 Hz. the computer 
solution included a uncorrelated noise source. The coherence shown in Fig. 7d is above the coherence 
threshold = 0.00639431 below 200Hz. Below 200 Hz. the coherence is above 0.15. Consequently, the 
cross spectrum phase angles look good below 200 Hz. Below 100 Hz., the uncorrelated noise level was to low 
to plot as shown in Fig. 7e. As shown in Fig. 7f, in order to match the measured phase angle the computer 
model puts the source at more than xa = 30 feet for frequencies less than 100 Hz. and less than xa = 30 
feet for frequencies greater than 100 Hz. 


NASA/TM— 2006-2 14352 


V. Discussion 


These results show how important it is to start with no preconceived idea of what acoustic source mode best 
represents the measured turbojet core noise. One should use a model that best fits the data and not ones 
idea of the model. The results show the uncorrelated noise source is not important at frequencies less than 
800 Hz at these angles. Furthermore, the results show that core noise dominates at frequencies less than 
800 Hz at these angles. Furthermore, at many operating points, the sound radiation field can be attributed 
to two coherent core noise sources. The results indicate that one might be the engine case while the other 
is the nozzle. The variation of the two-point model correlated source amplitude shown in Fig. 4d and the 
source separation results shown in Fig. 4a have the same frequency dependence that the combustion modes 
have in that they seem to be in 200 Hz bands. This suggests the correlated noise is due to the combustor 
source. The variation of the two-point model correlated source amplitude shown in Fig. 5d and the source 
separation results shown in Fig. 5a have the same frequency dependence that the combustion modes have 
in that they seem to be in 200 Hz bands. This again suggests the correlated noise is due to the combustor 
source. Note however, that the cross-spectrum phase angle measurements shown in Figs 4c and 5c have a lot 
of dispersion. The need to satisfy this dispersion may be causing the two correlated source with uncorrelated 
noise five parameter model to perform better than single correlated source with noise three parameter model. 
Consequently, further studies need to be done to resolve this issue. 

VI. Conclusions 

A new method to separate correlated signals buried in turbofan engine core noise has been presented. The 
method is based on finding acoustic model coefficients that enable a system of equations based on one or two 
correlated noise sources and an uncorrelated noise source to reproduce experimental data. The effectiveness 
and the reliability of the method have been tested using Pratt and Whitney PW4098 far field acoustic 
measurements. The agreement achieved between the experimental data and the acoustic model used confirm 
the effectiveness of the procedure. 
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Phase angle standard deviation, degrees 



Coherence 


Figure 3 . Standard deviation of phase angle of G34 based on 7| 4 and n 0 . 
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A. One-sided spectrum functions for acoustic signals from a point source 


We shall assume a compact source region and assume the source produces a wave which spreads spherically 
outward with no preferred direction. The wave equation in this case as discussed by Morse 30 page 309 is 

( 2 dp\ = ]_cPp . . 

r 2 dr \ dr ) c 2 dt 2 

We shall assume a simple source solution to the wave equation p(t) = exp(— ikR) sin(27r/ 0 t) where 


k 

cl 

3? 

MW 

TF 


2tt/ 0 


"ftftTK m 
MW ’ sec 

8.314, Joule 


0.02897 
80 


mole K 
kg 


mol 


(7) 

( 8 ) 


To get the auto-spectrum and cross-spectrum of a point source we use the following relationships from 
correlation and spectral analysis texts by BendatA 6 
The autocorrelation function of x is 


Rxx (t) = E[x k (t)x k (t + t)] 

The auto correlation function for a sum of two processes y(t) = x\ (t) + x 2 (t) is 


(9) 


Ryyir) = E[(x±(t) + + r) + x 2 (t + r))] 

Rx\X\{r^) T R Xi X 2 (t) ^X 2 X\ ("7") T Rx 2X2 (^") 

The two-sided auto and cross spectral density functions are 


Sxx(f) 


Sxy{f) 



Rxx (r) exp(-j2Trft)dt 
Rxy (t) exp(-j2nft)dt 


and of course we have the inverse relationships 


Rxx{r ) 


Rxyi'c) 


/ oo 

Sxxif) exp{+j2nfr)df 

-OO 

/ OO 

S xy (f)exp(+j2nfT)df 

-OO 


The two-sided auto spectrum function for a sum of two processes y(t) = X\ (t) + x?{t) is 


( 10 ) 


( 11 ) 


(12) 


Syy(f) = F[E[(xi(t) + X 2 (t))(x 1 (t + t) + X 2 (t + t))]] 

— E[R XlXl (t) -f- R XiX2 (t s ) T R X2X1 (t') -{- R X2X2 (t^ 

= S x1Xi (t) + S Xi x 2 (t) + S X2Xl (r) + S X2X2 (r)] (13) 

The auto-correlation function is the Fourier transform of the correlation of x(t) with itself is related to 
the Fourier transform of x(t) by 
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/ oo 

Rxx ( t ) exp(-j 27 r/r)dT 

-CO 

= X(f)X*(f) 


( 14 ) 

The cross-correlation function is the Fourier transform of the correlation of x(t) and y(t) is related to 
the Fourier transform of x(t) and y(t) by 


/ CO 

R xv {t) exp(-j 27 r/r)dr 

-OO 

= X(f)Y*(f) 


The one-sided spectral density function is G x (f) = 2 S x (f) 0 < / < oo . 

The auto-correlation function of a sine wave x(t) = A cos(27t /of) is 

A 2 

Rxx = y cos(2tt/o) 

The power spectral density of a sine wave x(t) = Acos(27r/ O t) is 

Gxx = y S(f - / 0 ) 

The one-sided spectrum function for a sum of two processes y(t) = X\ (t) + X2 (t) is 


(15) 


(16) 


(17) 


Gyy(f) — G XiXi (t) + G XiX2 (t) + G X2Xi (t) + G X2X2 (t)\ (18) 

Next, these relationships are applied to signals produced by acoustic waves from a point source and received 
at microphone i and j . 


Pat (t) 


P aj ( t ) 


exp(ikR ai ) sin( 27 r/oi) 

IX'di 

A 

-—exp{ikR aj )sm(2TTf 0 t) 


We have for the one-sided auto spectral density and the one-sided cross spectral density, 


GpatPai 
Gp ai Paj 




1 A 2 

2 R ai R a 


_ e -MR ai -Raj) §{f _ fo) 


B. Correlated noise source A 


Consider the i’th and j’th microphones and a single point acoustic source. 
As shown in Appendix A, the auto spectra are 


and the cross spectrum is 


Cp ai p ai — 2 ^2 /o) 

G Pa j Pa j = 

a j 


G 


Paj Paj 


1 A 2 

2 R ai Ra 


e -jk(R ai -Rt 




(19) 

( 20 ) 
(21) 


(22) 

(23) 

(24) 
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C. Correlated noise source A and uncorrelated source N 

Consider the i’th and j’th microphones, a single point acoustic source A and a noise source N. 
As shown in Appendix A, the auto spectra are 


G Pa .p a . 




and the cross spectrum is 


G 


Pa i Pa j 


1 A 2 

2 R ai Ra 


,-jk(R ai -R c 




(25) 

(26) 


(27) 


D. Correlated noise source A and correlated noise source B 

Consider the i’th and j’th microphones, a point acoustic source A and a point acoustic source B. 
As shown in Appendix A, the auto spectra are 


Gp a .pa. 

G P-jP-i 


1 A 2 
2^1“ 
1 A 2 


~ fo) + \§T S (f ~ fo) 


2Rl' 5 ^ ^ + 2 R 2 . 5 ^ ^ 


and the cross spectrum is 


G 


PiPj 


1 A 2 

2 Rdi Ra 


z -MR* t -R* j ) S ( f - fo ) + - 


-^- e - jk(Rbi - Rb ^S(f - fo) 

tibi ^bj 


(28) 

(29) 


(30) 


E. Correlated noise source A, correlated noise source B, and uncorrelated 

source N 

Consider the i’th and j’th microphones, a single point acoustic source A and a noise source N. 

As shown in Appendix A, the auto spectra are 


G n 


G 


Pa, Pa, 




1 A 2 

2 R 2 , 


- f 0 ) + - fo) + ~6{f - fo) 


1 N 2 

2 R 2 , 


and the cross spectrum is 


G 


PiPi 


1 A 2 

2 R a R a 


. e -MRa t -R* i ) S (f_ fo) + 


1 B 2 

2 Rb,Rb , 


^-jk(R bi -R b ) 


'*’S(f-fo) 


(31) 

(32) 


(33) 
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Figure 4. 582 rpm (N1 Corr.). 
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Cross spectrum Phase angle in degrees SPL re. 20 uPa. 
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Figure 5. 1622 rpm (N1 Corr.). 
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Figure 6. 1900 rpm.(Nl Corr.) 
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